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Abstract 

We take into account a generic form of a Dirac fermionic dark matter (DM) , which communicates 
with the Standard Model quarks via a scalar mediator, in a model-independent way. Four special 
interaction scenarios are investigated, where one is parity conserving and the other three are parity 
violating. Three of them result in the v suppressed DM-nucleon cross sections, where v ~ 10~'^c 
is the velocity of the DM in the laboratory frame. We constrain the masses of the dark matter 
and mediator as well as the couplings from the thermal relic abundance, and the recent results of 
the XENONIOO direct detection and collider experiments involving the two channels: (i) monojet 
plus large missing transverse energy, and (ii) dijet. The current monojet constraint is not stronger 
than that from the requirement of the correct relic density and the null result by the XENONIOO 
direct detection. We find that the dijet resonance measurements can exclude a large portion of 
the parameter space (m-^,my), where the couplings for the mediator coupled to the dark matter 
and to the quarks are small and have roughly the same magnitude. The constraint from indirect 
detections and diphoton resonance searches is also briefly discussed. 
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I. INTRODUCTION 



The nature of the dark matter is a challenging problem in the modern science. From 
astrophysical observations there are gravitational sources which are invisible to us besides 
the visible stars and galaxies. So far, we know little about the invisible objects and call 
them to be the dark matter (DM), which could be massive to gravitate and non-baryonic 
to keep the evidence of big-bang nucleosynthesis (BBN) intact. The DM is now believed to 
be responsible for ~ 23% energy density of our Universe, where the visible stars, galaxies 
and dim stellar gas only contribute ~ 4% [Ij. From the viewpoint of the particle physics, 
the DM may be a particle which goes beyond the standard model, and its identity remains 
elusive. One of the authors (K.C.Y.) has investigated the possibility that the galactic dark 
matter exists in an scenario where the phantom field is responsible for the dark energy [2] . 

The relic abundance can be used to determine the interaction strength between the ther- 
mal DM and the Standard Model (SM) particles. A stronger interaction keeps the dark 
matter staying in the thermal equilibrium longer, so that the Boltzmann factor further sup- 
presses its number density. To reproduce the correct dark matter relic [3], the governed 
Boltzmann's equation shows that the annihilation cross-section {av) of the dark matter into 
SM particles is of order 3 x 10~^®cm^/sec which is about 1 picobarn-c (pb-c) [4j. This 0{1) 
pb cross-section interacting with SM particles infers that the dark matter could be produced 
in recent collider experiments. The relic can be satisfied by the dark matter with mass being 
several hundred GeV, so that the dark matter can undergo a typical interaction with the 
electroweak scale cross section, which is the so-called weakly interacting massive particle 
(WIMP) miracle. Therefore it was suggested that the dark matter can couple to the Higgs 
boson in the effective theory below TeV scale and this kind of model is the so-called Higgs 
portal model [5]-[T2]. For the thermal dark matter, an upper bound on its mass has been 
pointed out in [13], where the DM mass < 34 TeV is set from the unitarity of the an- 
nihilation cross-section (av) with recent WMAP data. Some efforts are devoted to directly 
searching for WIMPs with masses of order < 10 GeV p3l [T5] . 

A number of underground experiments, e.g., XENON, CDMS and DAMA/LIBRA, have 
being performed to detect the DM directly scattered by the nuclei [T6HT9]. Although the 
controversial signals were detected by DAMA/LIBRA and CoGeNT groups, the null result 
has been reported by XENONIOO and CDMS groups, respectively. The XENONIOO data 
lead not only to the strongest limit so far for constraining the DM-nucleon cross section to 
below 10"*^^ cm^ for ~ 100 GeV (and 10""^^ cm^ for ~ 1000 GeV), but have also the 
constraints for cross section to below 10~^^ cm^ in the low mass region around 10 GeV [19]. 

The DM my be produced at the hadronic colliders. Since the DM interacts weakly with 
SM particles, they can escape from the detector. To detect the DM signals, additional SM 
particles are necessary for trigger in the production, so that a relevant process would be 
a combination of jets and missing energy. The process with monojet -|- missing transverse 
momentum final states has been reported by CDF [20], CMS [21], and ATLAS [221 123], 
and is one of the main channels for the dark matter searches at the hadronic colliders. 
Recently, ATLAS [22] has analyzed monojets with varying jet pt cuts using an integrated 
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luminosity of 1.00 fb ^. For all of these experiments, no obvious excess has been observed 
compared with SM backgrounds. 

On the other hand, if the dark matter interacts with the hadron via a mediator, it is 
possible to find out this mediator from the dijet mass spectrum at the hadron colliders, 
for which the searches by the CDF [21] and DO p5] collaborations have used data from pp 
collisions at the Tevatron, while searches by the ATLAS [26] and CMS [271 128] have used 
data from pp collisions at the LHC [21]. However, all measurements show no evidence for 
the new narrow resonance. 

In this paper, we take into account in a model-independent way that a Dirac fermionic 
DM can couple to the SM quarks via a scalar mediator, where the interaction can be parity 
conserving or violating. Adopting this framework, we will constrain the masses of the dark 
matter and mediator as well as the couplings from the thermal relic abundance, and the recent 
results of the XENONIOO direct detection and colliders involving the monojet measurements 
with large missing transverse energy [201423] and dijet resonance searches [21 ] [26H28] . Some 
works for the interaction between the DM and SM particles via a neutral spin-1 mediator 
can be found in Refs. [151 l30H33l 



In the effective Hamiltonian approach with a contact interaction between the DM and 
SM quarks, which is just suitable for the heavy mediator and has been discussed in the 
literature [33] - H^ . only an annihilation topology (usually the s-channel) contributes to the 
processes for the thermal relic abundance. However, for the general case, not only s-channel, 
but also the u- and t-channels may enter to participate the interactions, where the u- and 
t-channels can be switched off when the DM mass is smaller than the mediator or the s- 
channel is predominant if the coupling between the DM and mediator is much smaller than 
that between the quark and mediator. (See also Fig. [l]) 

The mass of the mediator for the DM interacting with SM particles could be comparable 
with the energy scale of the colliders, so that the interactions can be resolved and the 
mediator is produced on shell. In this condition, the description of the effective contact 
operator is no longer valid. There are two possible ways to detect the relevant effect at 
the colliders. One is to measure the monojet plus missing transverse energy in the final 
state, for which in addition to the SM background it is dominated by xx +monojet, with 
X being the dark matter. The other one is to search for the mediator Y in the dijet mass 
spectrum. We will show that the LHC monojet constraints on the DM mass as well as 
related parameters could be competitive with but not stronger than the dark matter direct 
detection at the present stage. For the dijet search, interchange with the mediator Y in 
the s-channel interaction is predominant in the bump like component of the resonances, for 
which the upper limit is set by the collider measurements. 

This paper is organized as follows. We begin in Sec. [IT] with a generic form of a La- 
grangian, which describes the Dirac fermonic DM interacting with the SM quarks via the 



scalar mediator. In Sec. HI, we perform a detailed numerical analysis for the relevant op- 
erators, and obtain the constraint on the parameter space of the DM and mediator masses. 
Together with discussions for the constraint from indirect detections and diphoton resonance 



searches, we conclude in Sec. IV 
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II. THE DIRAC FERMIONIC DARK MATTER MODEL 



In the present work, we are devoted to the study of the fermionic dark matter, which 
interacts with the standard model quarks via a scalar mediator (which might be exotic), in 
a model-independent way. The relevant parts in the Lagrangian are 

A£ = x{^?)- ^x)X + \{d,Yd^Y - m^F^) + YxiX^ + a^f)^ + Yf{\i + )/ , (1) 

where / is a unspecified SM quark, x the fermionic DM, and Y the scalar mediator which 
is the SM gauge singlet. The fermionic DM that we consider is assumed to be Dirac-like. 
(All the present calculations can equally apply to the Majorana case with a Z2 parity to 
ensure its stability.) The fermionic DM interacting with the SM particles through the scalar 
mediator is renormalizable and its origin is model-depedent. 

Our effective Lagrangian can be in connection with the Higgs portal model of the fermionic 
DM, for instance, for the following form of the potential [H llOr] 



A£ = -m^xx + YxiX^'s + *A^7')x + y + V2fH^H + -^{Y + V2){H^Hf + . . . , (3) 

where H is the SM Higgs doublet, for which the neutral component is shifted to {v + h)/ \/2 
with V = 246 GeV, and (f = Y + V2 is the real singlet scalar with the spontaneously breaking 
vacuum V2- The terms involving A4 and /i offer the Higgs portal between the dark and SM 
sectors. Comparing the Higgs portal model with our generic form, for the cases <C rrih, 
we have A{ ~ \iV2'fnf/m\ and A^ = 0. In the multi-Higgs doublet or next-to-minimal 
supersymmetric standard model (NMSSM) model |15], A^ can be non-zero, although it 
vanishes in the SM. 

For a long time, the model through Higgs portal is a nature generalization to connect the 
DM and SM sectors as studied in the literature. The Tevatron and recent LHC data hint 
that the SM Higgs may exist within the very narrow window nth ~ 125 GeV [IHHHOj. If it 
is true, the discovery of Higgs can claim the establishment of the SM of particle physics. 
Because the results of ATLAS and CMS are close to what we expect from the SM Higgs, it 
would exist little room for the invisible Higgs decays. In other words, the effective coupling 
of the Higgs to two F' s or to XX should be significantly small when the mass of the invisible 
final state particle is lighter than half of the Higgs mass. 

^ Another form used in the Uterature is [TTJ [37] 

A£ = -,7XoXX+%^H^i?XX. (2) 

For 171-^ <^ my, one can obtain Xh^-^/As — AJf2A4/my, and the mass of the DM x to be = mo — 
Ah^^z;^/(2As). An easier way to relate our resuh to the form of the Higgs portal model in Eq. ^ is to 
use the substitution in Eq. (jl|: r -> /i, A{ -> nif/v, A^ ^ 0, A^ = 0, and AJ -> X,,^^v/As. 
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III. NUMERICAL ANALYSES 



A. Dark matter constraint from the correct relic density 

The thermal rehc abundance of the dark matter particles after the freeze-out temperature 
can be solved to a good approximation [51] 



2 _ 1.07 X lO'^Xf ^ _ 1^ 0.038c/^mpim^(aann^rci) 

lLy.IL — , Xf — 111 , 

where h is the Hubble constant in units of 100 km/(s-Mpc), mp\ = 1.22 x 10^^ GeV is 
the Planck mass, is the dark matter mass, Xf = m^/Tf with Tf being the freezing 
temperature, g^: is the number of relativistic degrees of freedom with masses less than Tf, 
is the number of degrees of freedom of the x particle, and (cTann'i^rei) is the thermal average 
of the annihilation cross section, where Vj-e\ refers to the relative velocity of the dark matter 
particle during freeze-out. The current value for the DM density, coming from global fits of 
cosmological parameters to a variety of observations, is ^x^"^ ~ 0-112 ± 0.006 [Ij. 

The abundance of the DM is determined by the s-channel annihilation into a SM quark- 
pair through the exchange of the scalar Y, and, when E > my, by t- and w-channel anni- 
hilations into two Y particles with x niediated, where E is the energy of the dark matter 
particle. The processes are shown in Fig. [l} and the annihilation cross sections are listed in 
the following, 

^^^^^^^55 ^ ^/f pIp} , Emlplpy e{E - my) 



27r E^liAE"^ - mlY + mlTlsy] 27i [A{E^ - ml)ml + m^]^ 
.SP_^r ^"'^f P'xPf , Emlplpy e{E - my) 



271 E[{4:E^ - ml.)^ + mlTlpy] 2tx [A^E"^ - ml)ml + m\.f 
.PS_^r ^f^s' P) , mlplpl e{E - my) 



^""^'"'^ 27r E[(4E2 - m|-)2 + m^Vlsy] ^ 27r E'^[A{E^ - mf.)m2 + m^yf ' 

, ,PP _ Af Epf Xf mlplpl. 9{E - my) 

^ '''^ ' 27r [(4^2 - m2,)2 + m^Vlpy] ^ 2ti E^[A{E^ - m\)ml + mf.]2 ' 



where all flavors will be added on the right hand side, and p^iPv, and pf are the momenta 
of the X; ^ and quark, respectively. Here iVc = 3 is the number of the quark's colors. The 
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decay widths of the Y are given by 



with 



rss,Y = r(r ^ xx)s + r(r ^ ff)s , 
Tsp,Y = r(r -> xx)5 + r(r ^ /7)p , 
rp5,y = T{Y ^ xx)p + r(y ^ /7)s , 
rpp,y = T{Y -> xx)p + r(F ^ /7)p , (6) 



r(r ^ xx)p = ^myXfp^ e{my - 2m J , 
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- TV 9 

r(F ^ /7)P = -f myA^ , (7) 



where /3/ = ^1 — 4m^/my and = y''l — 4m^/my. Freeze-out happens at temperature 

~ m^/20, and thus the DM energy is around E ~ m^ + 3m^/40 (with Ty ~ m^w^/3). 
In Eq. (|5]), we have shown four special scenarios denoted by the superscript "5*5"', "SP", 
"P5", and "PP", for which only the couplings: (i) X} and A{, (ii) X} and A^, (iii) A^ and A{, 
and (iv) A^ and A^ are turned on, respectively. The scenario for SS is the parity conserving 
interactions, while the other three are parity violating interactions |^ The present relic can 
be approximated by 

ny- "iP^, (8) 

SO that it should be (dann'Wroi) ~ 1 pb ■ C. 

For the individual scenario given in Eq. (jsj), we shall consider three limits: (1) X}p = A{p 
for which the t- and w-channels are not negligible compared to the s-channel diagram, (2) 
20AJp = A{p for which the s-channel dominates over the t- and w-channels, and (3) AJ^ = 
100A{p for which the t- and u- channels predominate the contributions. To ensure that the 
calculation can be performed perturbatively, we impose the absolute value of A^ and A-^ 
to be smaller than 3. We shall consider the three conditions: (i) < Max[A-^, A-^] < 0.3, 
which is denoted by the magenta region, (ii) 0.3 < Max[A^, A-'^] < 1, which is given by the 
green region, and (iii) 1 < Max[A^, A-^] < 3, which corresponds to the blue region. We also 
assume the universal A{p, i.e, its value is independent of the flavor of the quark. It should be 
noted that for the (blue) region with A-^'-'^ > 1 the width of the scalar mediator Y becomes 
comparable with its mass, so that Y is no longer a good resonant state. 



If the mediator y is a pseudoscalar particle, then the PP scenario is parity conserving. 
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X fx Y X Y 

(a) (6) (c) 

FIG. 1. Tree level dark matter annihilation processes inferred from Lagrangian in Eq. 0. 



Under these scenarios, we show the allowed parameter space of (m^,my) in Fig. [2] 
Four remarks are in order. First, while maintaining the correct relic abundance, the region, 
where the contribution of the s-channel resonance is dominant, can be satisfied with small 
couplings, = A-^ and < A^^-^^ < 0.3. Second, the shaded region below the dashed line 
corresponding to E = my receives contributions from all channels, while above the dashed 
line, where + 3m^/20) < my, only the s-channel contributes and the t- and u- 

channels are not kinematically allowed. Third, for A^ ^ A-'^ < 1, the s-channel will be 
predominant; moreover, due to the nonzero width of the Y, the s-channel resonance region 
with a larger value of my is excluded, so that the allowed parameter space corresponds to 
the lighter mediator and lighter m^ is preferred in SS and SP scenarios. See also the central 
panels in Fig. |2} Forth, for A^ ^ A-^, the u- and t-channels are predominant in the region 
with my < m^, whereas the resonant region (with 2E ~ my), receiving only the s-channel 
contribution, is still allowed. See also the right panels in Fig. [2] 



B. Dark matter constraints from the direct detection and correct relic density 

The DM-nucleon interaction occurs via the exchange of scalar particle Y between the 
dark matter x and the nucleon in the t- channel process ■ Fo^^ the SS and PS 

scenarios, where A{ 7^ and A^ = 0, the interactions are spin-independent (SI) on the nucleus 
side, while for SP and PP scenarios, where A{ = and A^ 7^ 0, they are spin-dependent 
(SD). The elastic cross section can be expressed as 



^ef (xiV ^ xN) 
^fixN ^ xN) 



mzm% 



s s 
TT 


/ 


s p 


2 2 


7T 


2(m^ + mNyrriY^ 




2 2 
P mj^ 


TT 




A^'A^' 
'p 


4 

p 



vr 3(m^ + mjv)^my 



7 



ss , v=\i' 



SS . 20X^=1'' 



SS . A*=100Af 
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FIG. 2. Allowed parameter spaces of and my in reproducing the thermal abundance, (crvrei) ~ 
1 pb • c. The magenta, green, and blue regions correspond to < Max[A^,A^] < 0.3, 0.3 < 
Max[A^, A-'^] < 1, and 1 < Max[A-^, A-'^] < 3, respectively. From left to right, (i) A^^p = Xl^p for which 
the S-, t- and u-channels are equally important, (ii) 20As,p = A^^p for which the s-channel dominates 
over the t- and n-channels, and (iii) A^^p = lOOA^^p for which the t- and u-channels predominate 
the contributions. The shaded region below the dashed line receives contributions from s-, t- and 
ti-channels, while above the dashed line only the s-channel contributes. 
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where p = fiv with /i = m^mM / {my, + mAr) and v ~ 10~^c being the present velocity of 
the DM in the galactic halo. Here the effective coupling /tv of the nuclear matrix elements 
induced by the scalar SI coupling to quarks is given by 

/»= E — . (10) 

^-^ ^ nig 11 ° ^-^ mn 

q=u,d,s ^ Q=c,b,t ^ 

with f!j^^ and being hadronic matrix elements [52I455] . defined by {N\mqqq\N) = 

m^fi^fuNUN and'/g) = 1 - E,=u4,s fr? , where uj, is the spinor for the nucleon N and 
satisfied with the normalization unUn = irriN. We adopt the values used by the DarkSUSY 
package [55] . 

Ap) 



f^J = 0.023, /^^^ = 0.034, f^'J = 0.14, f^'> = 0.803 



4:^ = 0.019, 4:^ = 0.041, /g) = 0.14, fjff = 0.8. (11) 



The uncertainties of the hadronic matrix elements, which are mainly due to the magnitude 
of pion-nucleon cross section a^^j^f, might equivalently result in a factor two uncertainty in 
the product of the couplings, A^A-^. For obtaining the effective coupling of the nuclear 
matrix elements induced by the pseudoscalar SD coupling to quarks, we can perform the 
substitution 



{N'\qz-f,q\N) ^ (1 - r;5,)Ag(^)^«V75«7V , (12) 



where 6u = l,Sd = z, and 6s = w with r] = {1 + z + w)~^,z = niu/md, and w = mu/mg. 
Here Ag(^) measures the fraction of the spin carried by the quarks, q and g, in the nucleon, 
N. Thus gN is given by [56] 

g,, = {l- r/)A«(^)^ + (1 - r^^)Arf(^)^ + (1 - r/^)As(^)^ . (13) 
rriu rrid 

Again, we adopt the DarkSUSY numbers, 

Am(p) = Arf(") = 0.77, 
Ad^P) = An(") = -0.40 , 

As(P) = As(") = -0.12. (14) 

As Osi ~ (xx)(^Q') and Osd ~ (X7/i75X) (^7^750') operators dominate the interaction of 
DM with nuclear targets, the direct detection measurements quote the results as bounds 
on SI and SD cross sections per nucleon, respectively. However, it is interesting to note 
that numerically we have a^iixv ^ XP) — ^eiix''^ ~^ X'^) the SI and SD operators 
considered in this paper. Moreover, because gp and gn have the same sign for each case |^ 
these pseudoscalar SD interactions therefore behave as the SI one. 

^ Using the results given in |57J, gp and (?„ will have opposite signs, so that the resultant XENONIOO 
exclusion limit becomes even weaker in constraining the SP and PP scenarios. 
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In the direct detection, compared with erf; , one has that aff and o"^ are velocity sup- 
pressed by {v / cYw?^/ {m^ + m]yY and (f/c)^m^/(m^ + rriiy)'^, respectively, while cr^-^ is 
further suppressed by {v / cym%m'^/ {m^ + rriNY. As for ^ niN, we have a^f ^ a^/^ ^ 
a^i^ 3> (J^i^. These results are indicated in Fig. [s} 

Requiring that the correct relic density is obtained by the thermal freeze-out, and the 
DM-nucleus elastic cross section is consistent with the current 90%-C.L. upper bound by 
the XENONIOO measurement [11], we depict the allowed parameter space of (m^,my) in 
Fig. [sj For the "5*5"' scenario, in addition to the condition A-^ ^ for which the s-channel 
is predominant in the thermal relic density and few solutions exist, only the region that is 
very close to the s-channel resonance of annihilation in Eq. ([s]) is allowed. (See also Fig. |4j) 
If the Higgs-portal models of DM generate the operator as our "S'5'" scenario, the allowed 
solution will be very close to the s-channel resonance of annihilation, so that ~ 
[inilTT]. This is also called the "resonant Higgs portaT model [11]. The result will be further 
constrained by the invisible Higgs decay branching ratio Brj„i,, for which Brj„^ < 20% if 
the Higgs production rate at the LHC [IBtfKi] arises from the SM Higgs, or Brj„^, < 65% 
if new contributions exist [58]. As for other scenarios, compared with the constraint from 
the correct relic density, the direct detections offer marginal constraints on SP, PS and PP 
interactions. The PS scenario is related to the " pseudoscalar Higgs portaT model pTT] . 

Taking the thermal relic constraint and comparing the predictive curves of the dark 
matter-nucleon elastically scattering cross section to the current XENONIOO upper bound 
[19] . in Fig|4| one can easily realize that why the most portions of the (m^,my) region 
are excluded in the 5*5* scenario (except the resonance region) and in the SP scenario with 
A^ <^ A-^ , but the much greater parts are still viable in SP (with A^ = A-^ or A'^ ^ A-'^), 
PS and PP scenarios. For PS and PP scenarios, the predicted DM-nucleon cross section 
is fully below the XENONIOO bound and the main constraint comes from the thermal relic 
abundance. The results can be easily seen from Fig. |4]and are due to the fact that crff , cr^i^ 
and a^i^ are velocity suppressed, compared with af'^. Meanwhile, if the DM annihilation 
cross section is dominated by the t- and w-channel amplitudes, i.e., A'^ ^ A-^, the solutions 
mainly come from the region with > my, in addition to the narrow region close to the 
s-channel resonance where the couplings are small to maintain the correct relic density. 

Before proceeding, we remark that, as shown in Fig. [4], aei increases as becomes a 
much larger value for a fixed my and for A-'^ ~ A-^; meanwhile, if choosing a larger my as 
the input, it will result in a smaller aei- The reason is that, in the large m^ limit, because 
((Jann^^rei) ~ 1 pb ■ c (=const.) oc (A^A'^)^/m^, we therefore have the DM-nucleon scattering 
cross section to be aei oc (A^A''^)^/my oc mVmy. 
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FIG. 3. Same as Fig. [2] except for requiring that the correct rehc abundance is obtained and the 
DM-nucleon cross section aei is less than the XENONIOO upper Umit. The exclusion limit of the 
XNONIOO experiment has been extended to the region with > 1 TeV, assuming a linear fit in 
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FIG. 4. dark matter- nucleon elastically scattering cross section as a function of m^., predicted by 
requiring the values of couplings to be satisfied with the correct relic constraint, where we have set 
the bound of the couphngs to be Max[A'^, A-^] < 1. The XENONIOO experimental 90%-C.L. upper 
limit is depicted by the solid curve. From left to right, the long dashed, short dashed, dot-dashed, 
dotted curves respectively correspond to the following values, my = 100, 500, 1000, and 4000 GeV, 
as input. 
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C. Dark matter constraints from the LHC as well as the correct relic density and 
direct detection 

In addition to the DM mass constraints due to the correct rehc density and the direct 
detection, we will further consider the collider bound due to the monojet measurements with 
large missing transverse energy (-^r) in the final states, and dijet resonance search. 



1. Monojet + search 

At the hadron colliders, the dark matter particles can be produced and subsequently 
escape from the detector since they weakly interact with SM particles. Thus the signature 
of dark matter is a process with some missing energy. We can therefore use the monojet 
searches to constrain our dark matter model here. 

In experiments, monjet + final states have been studied by CDF [20], CMS [2T], 
and ATLAS [221 ES]. Recently, ATLAS [22] has analyzed monojets with varying jet pr cuts 
using 1.00 fb~^ of data, for which the number of the observed monjet events is in agreement 
with the Standard Model predictions. As studies in [311 EH] , the result of the very high pt 
(veryHighPt) selection region defined by the monojet study of ATLAS is the most effective in 
setting bounds for interactions involving the dark sector except for that the CDF data could 
provide a little more stringent bounds than LHC for < 25 GeV [60|. The ATLAS search 
in an event sample of the veryHighPt region has given 95% CL cross section upper limits of 
0.045 pb for new phenomena [22]. The leading channels containing the dark matters in the 
final state with the monojet plus missing transverse energy search at hadron colliders are 
shown in Fig. [sj while the Standard Model background consists mainly of [W — )■ £™^z/) + jet 
and {Z — 7- ujj)+iet final states, where denotes for the charged lepton lost by the detector. 
As the DM mass increases, the resulting production cross section decreases due to the smaller 
phase space, and therefore we can expect that the monojet constraint becomes weaker. 

Together with the constraints from the correct relic density and the exclusion bound of 
the XENONIOO direct detection measurement, we will use the recent ATLAS results on 
the monojet search to find the allowed parameter space of (my,m^). We have used the 
MadGraph 5 [61] and CTEQ6.1L parton distribution functions [62] to simulate the monojet 
plus missing energy events in the veryHighPt region, where the missing energy is > 300 
GeV, and one jet has a large transverse momentum pxiJi) > 350 GeV and pseudo- rapidity 
|r7(ji)| < 2. The detailed definition of the veryHighPt region can be found in [22] . 



2. Dijet resonance search 

The hadron colliders have performed the search for new particles beyond the standard 
model in the dijet mass spectrum. The process relevant to the dark matter search in the 
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FIG. 5. Dark sector productions through processes involving a single jet, where the jet can be one 
of quark, antiquark, and gluon, and the remaining two are from different accelerated proton beams 
at the LHC. 



experiments is the s-channel production and decay of the narrow dijet resonance, Y, shown 
in Fig. [6} where the contributions arising from the t- and u-channels are negUgible and will 
be discussed later. Note that the gluon-gluon fusion into Y through a quark triangle loop 
is highly suppressed, compared to qq annihilation into Y that we have shown in Fig. [6j 
CDF collaboration [21] using proton-antiproton collision data corresponding an integrated 
luminosity of 1.13 fb~^ has presented a search for new narrow particles whose decays produce 
dijets with invariant mass in the region 260 GeV < my < 1400 GeV. Recently, CMS and 
ATLAS collaborations [2S1 have published the dijet searches for new narrow resonances 
in the region 1 TeV < my < 4.3 TeV and 900 GeV < my < 4 TeV using integrated 
luminosities 5 fb~^ and 1 fb~^, respectively, at a center-mass-energy ^/s = 7 TeV. All the 
results show no evidence of new narrow resonance production over the SM background. 

For the dijet resonance search, the LHC data measured by the CMS and ATLAS collabo- 
rations give the most stringent constraint in the region 900 GeV < my < 4.3 TeV, while the 
CDF measurement gives stronger constraint for my < 900 GeV. Because the results given 
by CMS and ATLAS are quite similar, in the numerical analysis we will therefore mainly 
use the CMS and CDF data. The contributions arising from t-and u-channels in Fig. [6] can 
be negligible according to the following two reasons. First, the t- and w-channels mainly 
contribute to the processes for which the constituent (anit-)quarks are scattered with large 
transverse momenta, and can be highly suppressed by the pseudo-rapidity cut. Second, only 
the s-channel is relevant to bump like component of the resonances for which the upper limit 
is set by the collider experiments. See also the discussion in Ref. |29j . 



3. Results and discussions 

The numerical results are depicted in Figs. [7} [8] and [9j The monojet results are summa- 
rized as follows, (i) The monojet constraint is basically not stronger than the combination 
of the correct relic density and direct detection. As shown in Fig. [9| for the region satisfying 
1 < A-^(= A-^) < 3 and only for PS and PP scenarios, a very small fraction bounded by 
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FIG. 6. Leading dijet processes in the s-, t- and u-channels, where the initial and final states contain 
two quark jets and the mediator is the Y. (a) is the s-channel, (b), (c) and (e) are t-channels, and 
(d) and (f) are u-channels. Dijet experiments do not distinguish between the quark and antiquark. 



the black curve can be further excluded by the LHC monojet constraint. Note that, for the 
region with 1 < A^*^-^-* < 3, although the calculation is perturbatively convergent, the decay 
width of the Y particle may be comparable with its mass, so that Y is no longer a good 
resonant state, (ii) In Fig. [7| we show that, for Max[A^, A-^] < 1, the LHC monojet bound 
is always weaker than the constraint due to the correct relic density, (iii) Comparing our 
direct detection calculations in Eq. ^ with the results obtained from the operators defined 
in the effective theory framework, 

= , (15) 

where Fi^2 = 1 or 75, we have my/(A^A''^)^''^ = A. However, this relation is no longer correct 
in the calculations of the relic density and LHC monojet bound. For the relic density, it can 
be correct only when taking the limits: <^ my and Fy <^ my. In the LHC monojet 
search, because the LHC can not only resolve the interaction but also produce the on-shell 
mediator, so that my/(A^A/)^/^ depends on the value of my as shown in Fig. [rj 

For the dijet resonance search in the CDF and CMS experiments [211 [271 EH [63] , the best 
bin width is somewhat larger than the dijet invariant mass resolution cr, e.g. in the CMS 

1 /2 

measurements, a = 0.045 + 1.3/m^-j with rrijj being the dijet mass [HS]; the half-width of the 
resonance is significantly less than the experimental Gaussian resolution a. Conservatively 
speaking, these experimental results can be used to constrain the mass of the new particle 
with the narrow width Tpt ^ O.lm^j, where F^^ and mpt are the decay width and mass of 
the resonance. In Fig. [sl we show the observed upper limit at the 95% confidence level 
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on a X B X A for qq (or qq, qq) resonance from the inclusive analysis, compared to the 
predictions for processes with the s-, t-, and w-channel interactions via the scalar resonance, 
Y . We calculate the the results using CTEQ6.1L parton distribution functions [62]. Here a is 
the resonance production cross section, B is the branching fraction of the resonance decaying 
into the jet-jet final state, and A is the acceptance for the kinematic requirements. For CDF, 
the acceptance requires that the rapidity of the leading two jets satisfies \y\ < 1. For CMS, 
the acceptance corresponds to the pseudo-rapidity separation A?7 of the two leading jets 
satisfies [Ar^l < 1.3, and the two jets are also located in \rj\ < 2.5; A ^ 0.6 for isotropic 
decays independent of the resonance mass. Comparing the region plot of {my, m^) given in 
Fig. [9] with that in Fig. [3| we find that the dijet measurements offer strong constraint on 
the cases of = A^, where a large portion of the region corresponding to smaller couplings 
is further excluded. 



IV. DISCUSSIONS AND CONCLUSIONS 



Before we summarize, let us briefly discuss the constraints from indirect detections and 
diphoton resonance searches. Dark matter annihilation in the galactic halo can generate 
observable signals. The annihilation cross section, which can also read from Eq. (|5|, is 
proportional to for SS and SP scenarios, but is proportional to f° for PS and PP 
scenarios, where v ~ 10~^c is the velocity of the DM in the galactic halo. Compared to 
the thermal relic cross section (erf) ^ 3 x 10~^^ cm^ s~^ after rescaling the DM's velocity 
to be the same temperature, because the velocity of the DM is about 1/a/TOc at freeze- 
out temperature, the FERMI-LAT measurement from observations of Milky Way dwarf 
spheroidal galaxies [M] and BESS-Polar II [BS] data have no constraints on the SS and SP 
interactions due to their velocity suppression, but disfavor the dark matter mass < 30 
GeV for PS and PP interactions. However, one should note that the recent FERMI-LAT 
search for DM in gamma-ray lines and the inclusive photon spectrum |66j does not constrain 
all the interactions that we consider in this paper. 

The results for the diphoton resonance analyses at the collider experiments have been used 
to search for the Higgs [67] and Randall-Sundrum graviton [68]. As for the search for Higgs 
decaying into diphoton at the LHC and Tevatron, the three main Higgs production channels 
are the gluon-gluon fusion, vector boson fusion, and vector boson associate production. In 
our case, the Y particle is produced mainly from the qq annihilation. Requiring that the 
couplings for the DM into the quark pair are constrained by the correct relic abundance, 
we find that the cross section for pp into 77 via Y is at least 100 times smaller than that 
via the Higgs. Therefore, the contribution to the diphoton production due to the scalar Y 
particle can be negligible. In other words, the diphoton resonance experiments cannot offer 
the efficient constraint to the mass of the Y. 

In summary, we have studied the experimental constraints, which are due to measurements 
from the correct relic abundance, direct /indirect detections, and colliders, on a scenario that 
a Dirac fermionic DM interacts with SM quarks via a scalar mediator in a model-independent 
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way. We respectively consider four special scenarios denoted as SS, SP, PS and PP, where 
the former one is the parity conserving interaction, while the latter three are parity violating. 
Our present study can apply to the case of the pseudoscalar mediator, for which the PP 
scenario is the parity conserving interaction. For each interaction, we take three limits for 
couplings, X^p = A{p, A^^ ^ A{p and A^p ^ A{p, in the analysis. The main results are 
summarized as follows. 



Requiring A-^'-^ < 1, for which the resulting width of the scalar mediator Y is smaller 
than its mass, the allowed {m^,mY) parameter space that can maintain the correct 
relic abundance is very small for the coupling satisfied with AJ^ <C A{p or AJ^ ^ X^p, 
as shown in Fig. [2 For X}p ^ A{p, besides a very small portion of the parameter space 
corresponding to the s-channel resonance is allowed, the relic annihilation processes 
are dominated by the u- and t-channels, corresponding to the region satisfied with 
my < m^. 

In the direct detection, the SS and PS interactions, where A{ 7^ and A^ = 0, is spin- 
independent on the nucleus side, while the SP and PP interactions, where A{ = and 
A^ 7^ 0, is spin-dependent. For the latter, because the effective couplings gp ~ gn for 
the mediator coupled to the proton and neutron, these pseudoscalar SD interactions 
can be treated as the SI one. Compared with the direct detection cross sections a^f , 
we have shown that a^i^ and cr^"^ are velocity suppressed by {v/cYm'^/^m^ + rriNY 
and {v / cYm\ / {m^ + m]\f)^, respectively, where the DM velocity v ~ 10~^c, while 
a^i^ is further suppressed by (f /c)'^m^m^/(m^ + rriN)'^. As for ^ ttim, we have 



a^f ^ a^f ^ a^i^ 3> cr^^. These results are indicated in Figs. |3]and|4| Compared 
with the requirement by the correct relic abundance given in Fig. [2| most portions of 
the (m^,my) region can be further excluded in the SS scenario (except the resonance 
region) and in the SP scenario with A-^ ^ A-^ by the XENONIOO null measurement. 

The current monojet constraint is not stronger than that from the requirement of the 
correct relic density and the null result by the XENONIOO direct detection. Only a 
very small fraction located in the region satisfied with 1 < A^(= A-^) < 3 can be further 
excluded in PS and PP scenarios. However, although the calculation is perturbatively 
convergent in that region, the decay width of the Y particle may be comparable with 
its mass, so that Y is no longer a good resonant state. 



• We have used the results from the dijet resonance search in the CDF and CMS ex- 
periments to constrain the mass of the scalar mediator Y with the narrow width 
Ty ^ O.lmy. Fig. [9] has shown the viable parameter space of (DM's mass) and 
my (mediator's mass), which is satisfied with reproducing the correct relic abundance 
and constraints by the null result at the XENONIOO direct detection and the bounds 
due to the collider monojet and dijet searches. We find that the dijet measurements 
offer strong constraint on the case of A-^ ~ A'^, where a large portion of the region 
corresponding to smaller couplings is further excluded. 
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FIG. 7. my/(A^A^)i/2 
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i) The curves with a peak are satisfied with the cor- 
rect rehc density. From left to right, the red, orange, blue, and brown curves correspond to 
my = 100, 500, 1000, 4000 GeV, respectively, (ii) The region above dashed curves satisfies the AT- 
LAS null search for monojet + large final states with the constraint Max[A-^, A-^] < 1; from up 
to down, the brown, blue, orange, and red curves correspond to my = 4000, 1000, 500, 100 GeV, 
respectively. When the constraint is satisfied with Max[A^,A'^] = 1, the curves are denoted as 
dot-dashed (horizontal) lines, (iii) The XENONIOO experimental 90%-C.L. upper limit is depicted 
by the black curve, below which the region is excluded. For curves on the left hand side of the 
peak, with my ^ mx and narrow width limit of the Y, the value of my/(A^A-'^)^/^ is equivalent 
to A. 
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FIG. 8. The theoretical prediction for qq qq via the scalar resonance Y in colliders as a function 
of the mediator (resonance) mass, where the CDF (left panels) and CMS (right panels) 95% CL 
upper limits on a x B x A for the dijet production of the types gg, qq or qq^ are denoted by 
open boxes. By requiring the values of couplings to be satisfied with the correct relic density and 
narrow width Ty ^ O.lmy constraints, the predicted solid, dashed, dash-dotted, and dotted curves 
(from left to right) respectively correspond to = 100,300,500,700 GeV in the left panels, and 
= 500, 1000, 1500, 2000 GeV in the right panels. 
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FIG. 9. Same as Fig. [3] except that (i) the small region with white color but bounded by the black 
line is excluded by the ATLAS null search for monojet + large final states, and (ii) the white 
color region with smaller couplings is excluded by dijet measurements at the colliders. 
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